Direct evidence for a gapless Z^ spin liquid by frustrating Neel antiferromagnetism 
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By direct calculations of the spin gap in the frustrated Heisenberg model on the square lattice, with 
nearest- (Ji) and next-nearest-neighbor (J2) super- exchange couplings, we provide a solid evidence 
that the spin-liquid phase in the frustrated regime 0.45 < J2/J1 < 0.6 is gapless. Our numerical 
method is based on a variational wave function that is systematically improved by the application of 
few Lanczos steps and allows us to obtain reliable extrapolations in the thermodynamic limit. The 
peculiar nature of the non-magnetic state is unveiled by the existence of 5* = 1 gapless excitations 
at fc = (tTjO) and (0,7r). The magnetic transition can be described and interpreted by a variational 
state that is built from Abrikosov fermions having a Z2 gauge structure and four Dirac points in 
the spinon spectrum. 

PACS numbers: 



Introduction - During the "Valence-Bond-Solid era" 
most of the community working on highly-frustrated 
magnets believed that quantum spin liquids could not 
exist as true ground states of microscopic models and 
some kind of valence-bond order would have taken place 
in non-magnetic insulators (thus leading to trivial band 
insulators). Now, we are presently living the more ex- 
citing "Quantum-Spin-Liquid era" , where a plethora of 
different spin-liquid states are proposed as ground states 
of various magnetic systems, both theoretically and ex- 
perimentally-1 The turning point was marked by the dis- 
covery that stable gapped spin liquids may be found in 
effective low-energy Hamiltonians, which are based upon 
the so-called quantum dimer models^ or strong-coupling 
expansions.^ Since then, three main directions can be 
identified to study quantum spin liquids. The first one is 
the definition of ad hoc Hamiltonians that can be exactly 
solved to have a cartoon picture of the exotic proper- 
ties expected in generic systems (e.g., topological degen- 
eracy and fractional excitations) i^i^ The second one is 
the classification of different spin-liquid states according 
to hidden symmetries (i.e., beyond the Ginzburg-Landau 
description); examples may be given by the projective- 
symmetry group;^ tensor states^i or cohomology.^^-P Fi- 
nally, the third and more pragmatic one is to perform 
numerical simulations on Heisenberg or Hubbard models, 
in order to gain evidence that stable spin-liquid phases 
may indeed exist fibr- 
in this Letter, we take the latter point of view and 
investigate the J\—Ji spin- half Heisenberg model on the 
two-dimensional square lattice by systematically improv- 
ing accurate variational wave functions, to obtain a reli- 
able estimate of the exact ground state, along with few 
relevant low-energy states. This procedure allows us to 
extract the spin gap and show that a gapless spin-liquid 
phase exists in the highly frustrated regime. 
The Heisenberg J1—J2 model is defined by: 



H = Ji 



i^j) 



Ji • Oo 



J2 Y. ^- 



(1) 



where Si = {Sf,Sf,Sf) is the quantum spin operator 
on the site i; (...) and ((...)) indicate nearest- neighbor 
and next-nearest-neighbor sites. Here, we focus on the 
case where both super-exchange couplings are antiferro- 
magnetic and consider clusters with N = L x L sites and 
periodic boundary conditions. 

In the unfrustrated case with J2 = 0, it is well estab- 
lished that the ground state has Neel long-range order, 
with a staggered magnetization that is reduced from its 
classical value, i.e., M ~ 0.307^^1^ For large values of 
J2, the ground state shows again a collinear magnetic 
order with pitch vector Q — (7r,0) or (0,7r). The inter- 
mediate regime, around the strongest frustration point 
J2/J1 = 0.5, is the most debated one, since the com- 
bined effect of frustration and quantum fluctuations de- 
stroys antiferromagnetism and leads to a non-magnetic 
ground state. However, the nature of this quantum phase 
is still controversial. Since the pioneering worksji^r— 
it was clear that the problem was terribly complicated: 
many states can be constructed with very similar ener- 
gies but very different physical properties, e.g., having 
dimer or plaquette valence-bond order, or being totally 
disordered with short- or long-range resonating-valence 
bond fluctuations. This is mainly due to the fact that 
the non-magnetic region of the J1—J2 model is rela- 
tively small and several generalized susceptibilities may 
be quite large, ^'^ indicating that the ground state is on 
the verge of various instabilities. In this context, there is 
a convincing evidence that a third-nearest-neighbor cou- 
pling J3 may drive the system into a valence-bond solidj^ 

Recent density-matrix renormalization group (DMRG) 
results sparked the desire of understanding the phase 
diagram of the J1—J2 model, suggesting the existence 
of a true spin-liquid phaseJ^ In particular, by consid- 
ering cylindrical geometry, results for the singlet and 
triplet gaps provided some evidence for a fully gapped 
Z2 state in the region 0.41 < J2/J1 < 0.62, without lo- 
cal broken symmetry. Moreover, the calculation of the 
so-called topological entanglement entropy 7 was found 



to be consistent with the expected vahie of 7 = ln(2) for 
a gapped Z^ spin hquid. The most natural description 
of a fully gapped state is given in terms of the Schwinger 
boson representation of the spin operators. ^^ By per- 
forming a full optimization of the many-body wave func- 
tion on small sizes, we showed that this kind of bosonic 
ansatz may qualitatively reproduce some of the DMRG 
results j^ However, while in the weakly-frustrated regime 
the bosonic ansatz has magnetic order and excellent vari- 
ational energy, for 0.45 < J2I J\ ^ 0.6 a state constructed 
with Abrikosov fermions instead of Schwinger boson has 
better accuracyi^^ In fact, in a forerunner paperj^i three 
of us showed that, within this kind of fermionic rep- 
resentation, it is possible to have a particularly accu- 
rate description of the ground state in the strongly frus- 
trated regime. By using the language of the projective- 
symmetry group (PSG),— our variational wave function 
(dubbed Z2Azzl3 in Ref. 0) has a Zi gauge structure 
(implying gapped gauge excitations) but gapless spinon 
excitations with four Dirac points. A gapless spin liquid 
with topological degeneracy has been also suggested by 
using projected-entangled pair states on cylindrical ge- 
ometry.^^ These results suggest that a gapless spin liquid 
may be competitive with the one proposed by DMRG 
calculations. 

In this Letter, we present numerical calculations based 
upon a systematic improvement of the fermionic Z2Azzl3 
wave function that allow us to extract (i) the ground 
state energy, (ii) the energy of the lowest S ~ 2 state, 
and (iii) the energy of a state with S* = 1 and fc — (tt, 0) 
[or (0,7r)], so to extract the information about the exact 
spin gap. The state with S* = 1 and fc = (tt, 0) is par- 
ticularly interesting, since it is certainly gapped in the 
Neel phase and it is not expected to play any important 
role in a gapped non-magnetic regime (while it is one of 
the gapless modes in the coUinear magnetic phase that 
appears for large Ji values). On the contrary, this state 
is gapless in the Z2Azzl3 ansatz (see below for the de- 
tails). One of the main results of this work is to show 
that this S* = 1 excitation becomes indeed gapless in a 
region around Jij J\ = 0.5, and, therefore, a spin liquid 
with gapless triplet excitations both at fc = (vr, 0) and 
(tt, tt) represents the most natural candidate between the 
two magnetic phases characterizing the small and large 
J2/J1 regimes. 

Numerical method - The staring variational wave func- 
tions are defined through the mean-field Hamiltonian for 
the Abrikosov-fermion representation of the spin opera- 
torsi^ 

y-MF = YlU,jcl^Cj^^ + h.c. 

where for each bond {i,j) there are hopping (iij) and/or 
pairing (r/ij) terms; the mean- field Hamiltonian may also 
contain on-site terms (i.e., chemical potential and/or on- 
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FIG. 1: (Color on-line) Energies per site for the 5 = ground 
state (a), the S = 1 state with k = (7r,0) (b), and the 5 = 2 
with k = (0,0) (c) versus the variance for J2/J1 = 0.5. The 
results with p = 0, 1, and 2 are reported for L = 6, 8, and 10, 
and with only p = and p = 1 for L = 14 and L — 18. The 
variance extrapolated results are also shown. 



site pairing). Given any eigenstate |^mf) of the mean- 
field Hamiltonian, a physical state for the spin model can 
be obtained by a projection of it onto the subspace with 
one fermion per site: 



|*«)=^g|*mf), 



(3) 



where Vg = Y[ii''^i,t ~ "j.J-)^ i^ ^^^ Gutzwiller projec- 
tor, Tii^^ — c\ ^Ci^a being the local density. Depending on 
the symmetry of the mean-field ansatz (e.g., the pattern 
of the ii.j's and the ry^j's), the projected state may de- 
scribe different spin liquids, having for example C/(l) or 
Z2 gauge structure and gapped or gapless spinon spec- 
trum.^ 

Here, we will consider a state that is obtained by tak- 
ing a real pairing rjxy (with dxy symmetry) on top of 
the C/(l) state with nearest-neighbor hopping t and real 
pairing 773,2 _y2 (with d^^-y^ symmetry). The d^y term is 
crucial to break the U{\) gauge symmetry down to Z2. 
Restricting this coupling along the (±2, ±2) bonds im- 
plies commensurate Dirac points at A: = (±7r/2, ±7r/2) in 
the mean-field spectrum; with this choice, the optimal 
variational state is found by projecting the ground state 

oiUMF^ 

Besides the ansatz for the ground state, within this for- 
malism, it is straightforward to have simple representa- 
tions also for excited states. In this respect, it is useful to 
consider a particle-hole transformation for the down elec- 



trons on the mean-field Hamiltonian ([2]), i.e. 



"iA 



-^ Ci,i, 



such that the transformed Hamiltonian conserves the to- 
tal number of particles. Then, the ground state is ob- 
tained by filling the lowest N orbitals, with suitable 
boundary conditions (either periodic or anti-periodic) in 
order to have a unique mean-field ground state. Spin 
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FIG. 2: (Color on-line) Spin gap for the S = 1 excitation at 
k = (tt, 0) (left panel) and the S = 2 excitation with k — (0, 0) 
(right panel) for the 6x6 cluster. Results for p = 0, 1, and 
2 Lanczos steps are reported, together with the extrapolated 
and the exact ones. 



excitations can be obtained by creating the appropriate 
Bogoliubov quasi-particles (spinons) and possibly switch- 
ing boundary condition. By limiting to states that can 
be constructed from a single determinant, here we con- 
sider a S = 2 state with momentum k = (0,0) and a 
S = 1 state with k = (tt, 0) or (0,7r). By performing 
Monte Carlo calculations, we are able to compute sepa- 
rately the energies of these three states, so to assess the 
spin gap of the J1 — J2 model. 

In order to systematically improve the variational wave 
functions, we can apply a number p of Lanczos steps: 



l*p) = 1 



p 

?n — 1 






I*.), 



(4) 



where am are p additional variational parameters. 
Clearly, whenever jvl/^) is not orthogonal to the exact 
ground state, |^p) converges to it for large p. Unfortu- 
nately, on large sizes, only few steps can be efficiently 
afforded: here, we consider the case with p = 1 and 
p = 2 {p = Q corresponds to the original variational wave 
function) j^ Furthermore, an estimate of the exact en- 
ergy may be obtained by the variance extrapolation. In- 
deed, for a systematically convergent sequence of states 
\^p) with energy Ep and variance Up, it is easy to prove 



that E„ 



Ecx + const X (T^ 



where En 



-p, ....^.^ ^p - {^p\n\^p)/N 

and cr^ = ((*p|h2|*p) - (*p|7^|*p)2)/Ar are the energy 
and variance per site, respectively. Therefore, the ex- 
act energy E^^ may be extracted by fitting Ep vs Cp, for 
p = 0, 1, and 2. 

Results - A systematic analysis shows that the best 
possible ansatz for the variational wave function of the 
form (j31) has a non-vanishing d^y pairing in the whole 
regime 0.45 < J2/J1 ^ 0.6. Here, both the static 
structure factor S{q) and the dimer-dimer correlations 
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FIG. 3; (Golor on-line) The S = 2 spin gap as a function of 
the system size for the variational wave function and the Lanc- 
zos extrapolation for two values of J2. The thermodynamic 
extrapolations are consistent with a vanishing gap within the 
error bars, i.e., A2 = -0.04(5) and -0.07(7) for J2/J1 = 0.5 
and 0.55, respectively. The DMRG results on 2L x L cylinders 
(with open boundary conditions along x and periodic along y) 
for the 5* = 1 excitation are also shown.— Exact results (stars) 
of the S = 2 gap and the lowest S = 1 gap on the 6x6 cluster 
(with periodic boundary conditions) are reported. 



do not show any evidence for the occurrence of or- 
dered stateSfii^ in agreement with the DMRG results 
of Ref. m. 

Let us start by showing the accuracy of our method 
for the ground state and the two excitations: 5* = 1 at 
k = (tt, 0) and 5 = 2 at fc = (0, 0). In Fig. [H we report 
calculations for J2/J1 = 0.5 and different sizes of the clus- 
ter. For L = 6, where the exact results can be obtained 
by Lanczos diagonalizations, our extrapolations are ex- 
tremely accurate. Moreover, for the ground state, our 
best variational p = 2 state gives E/Ji = —0.503571(3), 
while Ecx/Ji = —0.50381; remarkably, the Lanczos step 
procedure remains effective even for larger sizes, the dif- 
ference between the best variational state with p = 2 and 
the extrapolated being very weakly size dependent (for 
L = 10, the p = 2 energy is E/Ji = -0.497549(2), while 
the extrapolated one is E/Ji = -0.49781(2)). The same 
applies also for excited states, see Fig. [T] The almost 
perfect alignment of the Lanczos steps, together with the 
impressive accuracy obtained up to relatively large clus- 
ters, clearly indicates that the exact ground state should 
be essentially described by the starting Z2 gapless state. 

In Fig. [2l we show the results for the 5 = 1 spin gap 
Ai at fc = (tt, 0) and the S — 2 spin gap A2 for the 
6x6 cluster, in comparison with the exact results. Re- 
markably, our approach based upon a spin-liquid wave 
function gives excellent accuracy on Ai in the whole re- 
gion 0.2 < J2/J1 < 0.55. A similar accuracy is also 
obtained for A2 in the strongly frustrated region (i.e., 
0.4 < J2/J1 < 0.55) even though this is not a sim- 
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FIG. 4: (Color on-line) The 5=1 with k — (tt, 0) spin gap as 
a function of the system size for the variational wave function 
and the Lanczos extrapolation for two values of J2 (upper 
panel). The behavior of the extrapolated gap as a function 
of J1I J\ is reported in the lower panel, the line is a guide to 
the eye. The Lanczos extrapolated gap as a function of L for 
different values of J 2 are also reported in the inset. 



pie excitation since it involves four spinons. Instead, at 
J2/J1 ~ 0.6 the accuracy deteriorates because a first- 
order transition to the collinear magnetic state takes 
place in the thermodynamic limiti ^^i^^'^^ in this region, a 
quasi-degeneracy of levels in the energy spectrum occurs, 
leading to a reduced overlap between the variational wave 
function and the lowest exact eigenstateJ^^'^^ 

Then, we consider larger cluster and perform a size 
scaling of the gaps, see Figs. [3] and H] for the S* = 2 and 
S = \ with k = (tt, 0), respectively. For L > 6, the ex- 
trapolations obtained with two (p = and 1) or three 
(p = 0, 1, and 2) points are perfectly consistent (i.e., the 
three points lie along a straight line, see Fig. [T]) There- 
fore, we perform the computationally demanding second 
Lanczos step only for relatively small clusters (up to the 
L = 10), while we limit to the first Lanczos step for large 
clusters (up to the L — 18). 

The S — 2 gap is reported for two values of the frus- 
trating ratio J2/J1, together with the S = \ gap ob- 
tained by DMRG calculations of Ref. [11. We find that 
the Lanczos step procedure clearly reduces the gap on 
any size. In contrast with the DMRG picture, we have 
a clear evidence that the spin gap closes when L — >■ 00 
for J2/J1 = 0.5 and 0.55. Indeed, the values that we 



obtain in the thermodynamic limit are both compatible 
with a vanishing gap, i.e., A2 = —0.04(5) and —0.07(7), 
see Fig. [3] We want to stress that our calculations are 
done on square clusters, having all the symmetries of the 
infinite lattice, and periodic boundary conditions, while 
DMRG calculations employed cylinders with 2Lx L sites 
and open boundary conditions along x. A possible expla- 
nation for having a finite gap within DMRG is that this 
method favors low-entangled states with finite gaps. On 
the contrary, our variational approach is more flexible, al- 
lowing both gapped and gapless states. At the pure p = 
variational level, the best wave function of the form ([3]) 
is found to be gapless, its energy being the lowest one 
among all states constructed from Schwinger bosons and 
Abrikosov fermions for 0.45 < J2/J1 ^ 0.6;^^ moreover, 
by applying few Lanczos steps, the finite-size gap lowers 
with no evidence for a finite value in the thermodynamic 
limit. 

Finally, the S — 1 gap with k — (tt, 0) has been com- 
puted for various values of J2 and cluster sizes, see Fig.|4l 
This gap is finite in the Neel phase for small J2/J1, where 
the only gapless S = 1 excitations have k = (0, 0) and 
k = {tt, tt). Indeed, this is what is found for J2/J1 < 0.48 
when the Lanczos extrapolation is considered, even if 
the starting variational wave function is gapless before 
Gutzwiller projection. Remarkably, in agreement with 
the theoretical picture of the Z2Azzl3 spin liquid, this 
gap vanishes for the two cases we investigated within the 
spin liquid region: J2/J1 = 0.5 and 0.55 (before the tran- 
sition to the collinear magnetic phase, which occurs for 
J2/J1 ^ 0.6). We expect that the 5 = 1 gap at k — (tt, 0) 
closes for J2 -^ J2 with a non-trivial exponent smaller 
than one, which is however not possible to estimate with 
our numerical results. Nevertheless, by performing a lin- 
ear fit of our data, we can obtain an upper bound of the 
Neel to spin liquid transition, which can be located at 
Ji = 0.48(2). 

Conclusions - In conclusion, by using a particularly ac- 
curate variational state and a procedure based upon the 
application of few Lanczos steps, we showed that it is pos- 
sible to extract important information on the spin gap of 
frustrated spin models. In particular, we provided a solid 
evidence that the spin-liquid phase of the J1—J2 model 
on the square lattice is gapless and may be very well de- 
scribed by using a Abrikosov-fermion mean field with a 
Z2 gauge structure and gapless spinons with four Dirac 
points at A: = (±7r/2, ±7r/2). The latter statement is fur- 
ther supported by the occurrence of a vanishing S = \ 
gap at the non-trivial momenta k — (tt, 0) and (0,7r). 
Our calculations give the first direct evidence for the ex- 
istence and the stability of highly-entangled gapless spin 
liquids in SU(2) spin models. 
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TABLE I: p=0 





J2/J1 = 0.4 


J2/J1 = 0.45 


J2/J1 = 0.5 


J2/J1 = 0.55 


L = 6 5 = 
5 = 1 
5 = 2 


-0.52715(1) 
-0.49582(1) 


-0.51364(1) 
-0.48557(1) 


-0.50117(1) 
-0.47541(1) 
-0.46265(2) 


-0.48992(1) 
-0.46538(2) 
-0.44974(2) 


L = 8 5 = 
5 = 1 
5 = 2 


-0.52302(1) 
-0.50835(1) 


-0.50930(1) 
-0.49635(1) 


-0.49656(1) 
-0.48453(1) 
-0.48095(1) 


-0.48487(1) 
-0.47299(1) 
-0.46806(1) 


L = 10 5 = 
5 = 1 
5 = 2 


-0.52188(1) 
-0.51362(1) 


-0.50811(1) 
-0.50080(1) 


-0.49521(1) 
-0.48830(1) 
-0.48722(1) 


-0.48335(1) 
-0.47625(1) 
-0.47443(1) 


L = 14 5 = 
5 = 1 
5 = 2 


-0.52124(1) 
-0.51772(1) 


-0.50745(1) 
-0.50425(1) 


-0.49447(1) 
-0.49131(1) 
-0.49141(1) 


-0.48242(1) 
-0.47904(1) 
-0.47880(1) 


L = 18 5 = 
5 = 1 

5 = 2 


-0.52107(1) 
-0.51921(1) 


-0.50728(1) 
-0.50553(1) 


-0.49426(1) 
-0.49249(1) 
-0.49274(1) 


-0.48215(1) 
-0.48018(1) 
-0.48026(1) 





T 


ABLE IL p=] 


L 






J2/J1 = 0.4 


J2/J1 = 0.45 


J2/J1 = 0.5 


J2/J1 = 0.55 


L = 6 5 = 
5 = 1 
5 = 2 


-0.52928(1) 
-0.50042(1) 


-0.51538(1) 
-0.49020(1) 


-0.50323(1) 
-0.48082(1) 
-0.46807(1) 


-0.49303(1) 
-0.47238(1) 
-0.45605(1) 


L = 8 5 = 
5 = 1 
5 = 2 


-0.52501(1) 
-0.51157(1) 


-0.51101(1) 
-0.49963(1) 


-0.49855(1) 
-0.48857(1) 
-0.48489(1) 


-0.48777(1) 
-0.47847(1) 
-0.47305(1) 


L = 10 5 = 
5 = 1 
5 = 2 


-0.52368(1) 
-0.51610(1) 


-0.50973(1) 
-0.50344(1) 


-0.49718(1) 
-0.49165(1) 
-0.49041(1) 


-0.48622(1) 
-0.48090(1) 
-0.47867(1) 


L = 14 5 = 
5 = 1 
5 = 2 


-0.52287(1) 
-0.51966(1) 


-0.50899(1) 
-0.50632(1) 


-0.49638(1) 
-0.49398(1) 
-0.49387(1) 


-0.48519(1) 
-0.48270(1) 
-0.48221(1) 


L = 18 5 = 
5 = 1 

5 = 2 


-0.52259(1) 
-0.52083(5) 


-0.50874(1) 
-0.50137(1) 


-0.49611(1) 
-0.49475(1) 
-0.49485(1) 


-0.48475(1) 
-0.48327(1) 
-0.48319(1) 



TABLE III: p=2 





J2/J1 = 0.4 


J2/J1 = 0.45 


J2/J1 = 0.5 


J2/J1 = 0.55 


L = 6 5 = 
5 = 1 
5 = 2 


-0.52957(1) 
-0.50130(1) 


-0.51558(1) 
-0.49108(1) 


-0.50357(1) 
-0.48197(1) 
-0.46929(1) 


-0.49399(1) 
-0.47419(1) 
-0.45750(1) 


L = 8 5 = 
5 = 1 
5 = 2 


-0.52539(1) 
-0.51224(2) 


-0.51125(1) 
-0.50033(1) 


-0.49886(1) 
-0.48952(1) 
-0.48583(4) 


-0.48841(2) 
-0.48008(4) 
-0.47443(2) 


L = 10 5 = 
5 = 1 
5 = 2 


-0.5240(1) 
-0.51671(7) 


-0.51001(1) 
-0.50398(1) 


-0.49755(1) 
-0.49243(1) 
-0.49121(3) 


-0.48693(3) 
-0.4825(2) 
-0.4800(2) 


L = 14 5 = 
5 = 1 
5 = 2 










L = 18 5 = 
5 = 1 
5 = 2 











TABLE IV: extrapolation 





J2/J1 = 0.4 


J2/J1 = 0.45 


J2/J1 = 0.5 


J2/J1 = 0.55 


L = 6 S = 
S = l 
S = 2 


-0.52972(1) 
-0.50204(5) 


-0.51566(1) 
-0.49187(4) 


-0.50382(1) 
-0.48312(6) 
-0.4706(1) 


-0.49521(7) 
-0.4766(1) 
-0.4587(1) 


L = 8 S = 
5 = 1 
S = 2 


-0.52556(1) 
-0.51282(1) 


-0.51140(1) 
-0.50085(1) 


-0.49906(1) 
-0.49039(2) 
-0.48677(1) 


-0.48894(3) 
-0.48194(3) 
-0.47602(3) 


L = 10 S = 
S = l 

S = 2 


-0.52429(2) 
-0.51718(3) 


-0.51017(2) 
-0.50445(3) 


-0.49781(2) 
-0.49329(5) 
-0.49203(5) 


-0.48766(6) 
-0.4842(1) 
-0.48157(8) 


S = l 
S = 2 


-0.52351(2) 
-0.52052(2) 


-0.50953(1) 
-0.50724(3) 


-0.49722(2) 
-0.49562(5) 
-0.49539(4) 


-0.48696(5) 
-0.48594(7) 
-0.48524(9) 


L = 1S S = 
S = l 
S = 2 


-0.52333(1) 
-0.52180(4) 


-0.50940(1) 
-0.50828(3) 


-0.49717(2) 
-0.49645(3) 
-0.49636(3) 


-0.48698(5) 
-0.48656(5) 
-0.48638(5) 



